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This work examines the heat transfer analysis of a convective flow over a vertical plate under the combined influence of viscous 
dissipation and thermal radiation in the presence of heat source/sink with the plate being subjected to a variable surface 
temperature. The governing boundary layer equations are formulated, simplified and non-dimensionalised. The dimensionless 
equations are solved by employing Crank Nicolson’s implicit finite difference scheme. The effects of dimensionless numbers 
affecting the flow are shown graphically on the dimensionless velocity and temperature profiles. 


INTRODUCTION 

Free convective flow of a viscous, incompressible fluid past semi- 
infinite or infinite flat plates with thermal radiation effects is of great 
importance due to their numerous applications in engineering and other 
processes involving high temperatures such as nuclear power plants, gas 
turbines and thermal energy storage. It is often assumed that in such 
flows, the viscous dissipative heat is negligible. It was shown by 
Gebhart (1962) that the viscous dissipative heat is important when the 
natural convection flow field is of extreme size or the flow is at 
extremely low temperature or in high gravity field. Gebhart and 
Mollendorf (1969) also investigated the effects of viscous dissipative 
heat on free convective flow past semi-infinite plates. Analysis on the 
radiation effects of free convection flow of a gas past a semi-infinite flat 
plate was considered by Soundalgekar et al. (1960). Sakiadis (1961) 
studied the steady flow on a moving continuous flat surface and 
developed a numerical solution using a similarity transformation. 
Sparrow and Cess (1961) investigated the effects of magnetic field on 
the natural convection heat transfer. Romig (1964) considered the effect 
of electric and magnetic fields on the heat transfer to electrically 
conducting fluids. The heat transfer characteristics in the laminar 
boundary layer of a viscous fluid over a stretching sheet with viscous 
dissipation or frictional heating and internal heat generation was 
analysed by Vajravelu and Hadjinicolaou (1993). Muthucumaraswamy 
and Ganesan (2001) worked on the effect of the chemical reaction and 
injection on flow characteristics in an unsteady upward motion of an 
isothermal plate. The radiation and mass transfer effects on an unsteady 
MHD free convection flow past a heated vertical porous plate with 
viscous dissipation was examined by Gnaneshwara and Bhaskar (2009). 
Alam, Rahman and Sattar (2009), tried transient magnetohydrodynamic 
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free convective heat and mass transfer flow with thermophoresis past a 
radiative inclined permeable plate in the presence of variable chemical 
reaction and temperature dependent viscosity. Anjali Devi and Ganga 
(2010) considered the problem of dissipation effects on MHD nonlinear 
flow and heat transfer past a Porous surface with prescribed heat flux. 
(2009) investigated the radiative heat transfer with 
hydromagnetic flow and viscous dissipation over a stretching surface in 
the presence of variable heat flux. Gribben (1965) investigated the 


Kumar 


boundary layer flow over a semi-infinite plate with an aligned magnetic 
field in the presence of pressure gradient where he obtained solutions for 
large and small magnetic Prandtl number using the method of matched 
asymptotic expansion. Unsteady two-dimensional laminar free 
convection flow of an incompressible, electrically conducting fluid 
through a porous medium bounded by infinite vertical plane surface of 
constant temperature was carried out by Helmy (1998). The radiation 
effects on mixed convection along a vertical plate with uniform surface 
temperature using Keller Box finite difference method was investigated 
by Takhar (1996). Chamkha et al. (2001) examined the radiation effects 
on free convection flow past a semi-infinite vertical plate with mass 
transfer. The radiative and free convections effects on the oscillatory 
flow past a vertical plate was carried out by Mansour (1990). Raptis and 
Perdikis (1999) studied the effects of thermal radiation and free 
convective flow past moving plate. Gregantopoulos et al. (1981) 
considered two-dimensional unsteady free convection and mass transfer 
flow of an incompressible viscous dissipative and _ electrically 
conducting fluid past an infinite vertical porous plate. However, 
convective flow over a vertical plate under the combined influence of 
viscous dissipation and thermal radiation in the presence of heat 
source/sink with the plate being subjected to a variable surface 
temperature has not been considered so far in literature. 
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MATHEMATICAL ANALYSIS 


Consider an unsteady flow of an electrically conducting, viscous, 


incompressible and radiating fluid flowing past a semi-infinite vertical 
plate. The x axis is taken along the plate in the vertical direction while 
the y axis is taken normal to it. A uniform magnetic field with strength 
B, is applied perpendicular to the plate along the y axis and the effect of 
viscous dissipation is taken into account. All the fluid properties are 
taken to be constant except the influence of the density variation in the 
body force term that are caused by changes in temperature which are 
approximated by Boussinesq. The governing boundary layer equations 
are as follows 


du Ov 

a ae 0 (1) 

du du Ou * = 78 y 

se ta ee, shop = 1) (2) 
oT* oT* aT* a? T* 1 Oqr , Qs(T*-TS) m (2) 
at ae th Gy ~ * oy ne oy pCp seer ay 
(3) 


The initial and boundary conditions are given as 


u=0,v=0,T*=TZ forallx,y,t <0 
uU=Uy,v=0,T* =Ty, aty=O0fort>0 
u>0,T*>-TZ asy>o (4) 


where u is the velocity in x direction, v is the velocity in the y 
direction,t* is the time, g is the acceleration due to gravity,T* is the fluid 
temperature, Tf is the free stream temperature, PB is the coefficient of 
thermal expansion,ug is the velocity of the plate,v is the kinematic 
viscosity, By is the magnetic field, q, is the radiative heat flux, o is the 
electrical conductivity, C, is the specific heat at constant pressure, p is 
the fluid density, q, is the radiative heat flux and Q, is the heat 
source/sink parameter. 


For an optically thin limit, the radiating gas is said to be non-gray 
near equilibrium and it is defined according to Reddy et al. (2012) as 
follow 


eae = 4(T* —TA)L (5) 


d . : beat 
a), da, kA, is the absorption coefficient and ep, 


where L= if kay ( 


is the plank function. 


By introducing these non-dimensional quantities, 


x u v t*u, T*-TE te 
X=*,ys4,ust vel c= r= g= = 
l l Uo Uo l T= Teo PCpUo 
(6) 
_ oBRl = us _v _ ALL _ gBT Ty-Te)l 
M= jHe= eet! et Gea ee CO 
pu Cp(Tw-T. 2) a PCyUo Uo 


the dimensional equations (1), (2) and (3) becomes 


=0 (7) 
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dU aU dU 1 0°U 

ap oe ae ne ave + Gr? — MU (8) 

aT aT aT 1 @T au? 

at ag + ay some oye + Be (7) —NT—1) + QT 


(9) 
together with the dimensionless initial and boundary conditions 


U=0,V=0,T=0 forallx,Y,t<0 
U=1,V=0,T=1atY=Ofort>0 
U-0,T-0 


asY > 0 (10) 

Where U and V are the dimensionless velocities in X and Y directions, t 
is the time, T is the temperature of the surrounding fluid, Re is the 
Reynold number, G,. is the thermal Grashof number, M is the Magnetic 
field parameter, Q is the heat absorption/generation, Ec is the Eckert 
number, N is the thermal radiation parameter and Pr is the Prandtl 
number. 


Numerical Approach 

The set of unsteady non-linear partial differential equations under the 
boundary conditions are solved by employing an implicit finite 
difference scheme of Crank-Nicolson type. This method converges 
faster and it is unconditionally stable. The corresponding finite 


difference equations are given in the form below: 


uift— uy +uf K UE atu ue tue a -uk 12 “i ¥ vith Vike + +5 -Vi_ 1 
4AX 2AY 
0 (11) 
k+ k+ k+ k k+ K+1 
Ui; . “Hh 4 Uk, (Ug; *u7e Lit Uig— UE, ;) +V% (Ug ia- Uj, 4 +R gy uk; a) _ 
At 2AX 4AY 
1 Uk -2uit+uK tuk. —-Ui+U Kia Ti 47K ukt*+uh, 
+Gr —M (12) 
PrRe 2(AY)2 2 
k+1 k k+1 k+1 k k 
i Hoy a Shea a 
At bd 2AX 
k+1 k+1 k k 
+ ve (Tht ~ T;j 1+ Tij41 — Ta) = 
Tn 
k+1 k+1 k+1 k k k 
1 Te 20g Pigs + Tiga — Ty + pa 
PrRe 2(AY)2 
k k 2 
eee Uij4r — Vij 
2AY 


Titer, petty 
—N (OS ar) 1| +Q-—— TT (13) 
In the finite difference equations, the co-ordinates of the mesh points 
iAX,Y = jAY,t = kAt where 1, j and 
k are positive integers and the value of U and V at the mesh point is 


of the solution are defined by X = 


given asUjs, and Vi‘. 

The region of integration is considered to be a rectangle with sides 
Xmax=! and Ynqx=l whereby regarding Y;,.q,=1 as Y=oo which lie very 
well outside the momentum, thermal and concentration boundary layers. 
The subscript i, j and k relates to grid points along X, Y and t 
respectively. We therefore divide X and Y into M and N grid spacing 
respectively. The mesh sizes are AX=0.05, AY=0.05 and At=0.01. 

The coefficients Ui; and Viki appearing in the difference equation 


are treated as constants in any one time step. The values of T, U and V 
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Figure 1: Effect of Variable Temperature (Q) on Temperature Profile 
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Figure 2: Effect of Eckert number (Ec) on Temperature Profile 
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Figure 3: Effect of Prandtl number (Pr) on Temperature Profile 
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Figure 4: Efeect of thermal radiation (N) on Temperature Profile 
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are known at all grid point from the initial conditions at t=0. The 


computations off, U and V at time level (k+1) using the known values 
at previous time level k are carried out as follows. 

At every internal nodal point, equation (13) on a particular i-level 
forms a tri-diagonal system of equations which is solved by Thomas 
algorithm as discussed in Carnahan et al (1969). 

In this way, the values of T are found at every nodal point for a 
particular i at (k + 1) time level. Using T values at (k + 1) time level in 
equation (12), the values of U at (k + 1) time level are found in the same 
way we calculated T values. Thus, the values of T and U are known on a 
particular i-level at every nodal point. The values of V are calculated 
explicitly using equation (11) at every nodal point on a particular i-level 
at (k + 1) time level. This process is repeated for various i-levels. Thus, 
the values of T, U, and V are known at all grid points at (k + 1) time 
level. Computations are carried out until the steady state is reached. The 
steady-state solution is assumed to have been reached, when the absolute 
difference between the values of U as well as temperature T at two 
consecutive time steps are infinitely small. 


DISCUSSION OF RESULTS 

To report on this work, calculations are carried out for different values 
of fluid parameters to show their influence on the flow field. Default 
values are Gr=5.0, Re=1.0, M=1.0, Pr=0.7, N=0.04, Q=1.0, Ec=0.002. 
As observed from figure 1, the temperature distribution decreases as 
heat absorption parameterQ increases because when heat is absorbed, 
the buoyancy force decreases the temperature profiles. A similar 
observation is seen by considering various values of Eckert number Ec 
as shown in Figure 2. To observe the effect of Eckert number, other 
parameters are kept fixed, i.e Gr=5.0, Re=1.0, M=1.0, Pr=0.7, N=0.04, 
Q=1.0, n=0.5 and the effect of Eckert number is observed for Ec= 0.0, 
0.1, 0.4, 1. It’s graphical representation for the temperature distribution 
revealed that by increasing Eckert number Ec, the temperature profile 
increases. Figure 3 illustrates the numerical results for various values of 
Prandlt number Pr=0.04, 0.7,3. It is observed that by increasing Prandtl 
number Pr, the temperature distribution decreases. This is as a result of 
the presence of Lorentz force that has power to slow down the motion of 
the fluid. Figure 4 represents the influence of radiation on the 
temperature profile. Increasing radiation value reduces temperature 
profile. 
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